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In this article we use 1420 MHz data to demonstrate the likely reality of Galactic radio Loops V and VI. We further 
estimate distances and spectral indices for both these and the four main radio loops. In the cases of Loops I - IV, radio 
spectral indices are calculated from the mean brightnesses at 1420 and 820/404 MHz. The spectral indices of Loops V 
and VI are obtained from T — T plots between 1420 and 408 MHz. Using the supernova remnant (SNR) hypothesis for 
the origin of radio loops, distances are calculated from the surface brightnesses and the angular diameters at 1420 MHz. 
We also study how results for brightnesses and distances of radio loops agree with current theories of SNR evolution. For 
this purpose, the ambient density and initial explosion energy of the loops are discussed. We also discuss applications of 
different E — D relations. The results obtained confirm a non-thermal origin and nearby locations for the Galactic radio 
loops. Therefore, we have indications that they are very old SNRs that evolve in low ambient densities, with high initial 
explosion energies. 
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1 Introduction 

It is well known that some spurs seen in the radio contin- 
uum background distribution can be associated such as to 
trace small circles. The set of spurs which form a particu- 
lar small circle is called a loop. Four major loops had been 
recognized by the early seventies, these being discovered 
and studied in the following order: Loop I (Large, Quigley 
& Haslam 1962; Haslam, Large & Quigley 1964; Large, 
Quigley & Haslam 1966; Salter 1970), Loop II (Large et 
al. 1962; Quigley & Haslam 1965; Salter 1970), Loop III 
(Quigley & Haslam 1965, Salter 1970) and Loop IV (Large 
et al. 1966; Salter 1970; Reich & Steffen 1981). The most 
precise determination of parameters for these circles was 
made by Salter (1970) and published by Berkhuijsen, Has- 
lam & Salter (1971). Salter used the then best data avail- 
able at 408, 404, 240 and 178 MHz. A detailed review of 
the subject was published by Salter (1983). It was noticed 
by Milogradov-Turin (1970, 1972, 1982) and Milogradov- 
Turin & Urosevic (1997) that some other spurs could be 
connected into loops. They proposed Loop V, formed by 
negative latitude spurs in Taurus, Pisces and Pegasus, and 
Loop VI, corresponding to weak positive latitude spurs in 
Leo and Cancer. Loop V has well defined "end" spurs and 
weak "bridges" between them, while Loop VI is practically 
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equivalent to the high brightness temperature region betwe- 
en cold holes near the North Galactic Pole. They further 
computed parameters for the proposed loops using the sur- 
vey at 38 MHz by Milogradov-Turin & Smith (1973). 

The radio surveys at 1420 MHz (Reich & Reich 1986) 
and at 408 MHz (Haslam et al. 1982) are used as the basic 
source of data in the present paper. The northern-sky map of 
the radio continuum brightness at 1420 MHz (data are taken 
from Reich & Reich 1986) is presented in Figure [T] The 
paths of Loops V and VI are marked in order to show the 
approximate positions and dimensions of these loops. The 
small circles represent data from Table 1 of Milogradov- 
Turin & Urosevic (1997) which contains approximate pa- 
rameters for the best-fit circles. Although other theories ha- 
ve been suggested for the origins of the radio loops, (for 
example, Sofue 1994, 2000 suggested that the radio loops 
and the diffuse X-ray emission around the Galactic center 
can be explained by a giant explosion in the Galactic center 
region), we present results which suggest that radio Loops 
V and VI may represent supernova remnants (SNRs). 

Recent S — (surface brightness - diameter) relations 
(Urosevic 2002; Arbutina et al. 2004) have been used in this 
paper to estimate distances for the Galactic radio loops. 

2 The Data 

The survey at 1420 MHz (Reich & Reich 1986) provided 
a basis for the recomputing of small circles for the major 
loops and Loops V and VI. The observational data were ob- 
tained with the Stockert 25 m telescope. The effective sen- 
sitivity is about 50 mK in Tb (where Tb is for an average 
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Fig. 1 The 1420 MHz northern-sky continuum survey. 
The contour lines represent brightness temperatures of 
3400, 3600 and 3800 mK. The spurs in Taurus, Pisces and 
Pegasus form Loop V and the spurs in Leo and Cancer form 
Loop VI. The loop centers are marked by small white cir- 
cles. The temperature scale in mK is given below the figure. 

brightness temperature), and the angular resolution is 35'. 
The survey at 408 MHz (Haslam et al. 1982) is an all-sky 
survey obtained by combining surveys from four telescopes, 
the Jodrell Bank Marks 1 and lA, Effelsberg and Parkes, 
with effective sensitivity of about 1.0 K in Tb and an an- 
gular resolution of 51'. These data are available from the 
MPIfR Survey Sampler ("Max-Planck-Institut fiir Radio- 
astronomie", Bonn), which is available on-line at: http:// 
www.mpifr-bonn.mpg.de/survey.html and allows users to 
pick a region of the sky and obtain images at a number of 
wavelengths. 

3 Analysis 

One way of understanding SNR evolution is by studying 
the relation between the surface brightnesses and diameters 
of SNRs, the so-called E - £) relation. Shklovskii (1960) 
developed the first theoretical Y, — D relation, which is as- 
sumed to be of the form: 

E = ^D-'^, (1) 
where D is the diameter of the object, and A is a coeffi- 
cient which does not depend on D, but does depend on the 
initial energy that the supernova provided to the remnant 
and on the density of unperturbed matter through which the 
remnant is expanding. 7 is a function of a, the flux density 
spectral index [Si, oc i/^"). 

The first empirical Yj — D relation was derived by Milne 
(1970), while Berkhuijsen (1973) was the first person to in- 
clude the Origem Loop and Loop I in the S — relation. 
All the main radio loops (Loops I, II, III and IV) were in- 
cluded in the study of the Y — D relation of SNR at known 
distances by Berkhuijsen (1986). 



Brightness temperatures over the areas of the loops, as 
required for our T — T plots, were obtained from the con- 
tinuum survey data of Reich & Reich (1986; at 1420 MHz) 
and Haslam et al. (1982; at 408 MHz). Each loop area was 
divided into different sections (corresponding to spurs) and 
estimates for these sections were combined. The areas used 
for the individual loops were not defined with respect to 
their best-fit circles, but were obtained from the radio con- 
tinuum maps. Assuming that Loops I- VI are old SNRs, pro- 
bably in their isothermal (radiative) phase, only the data 
from the relevant spurs of the loops were used for deriv- 
ing brightness temperatures. The surface brightnesses of the 
SNRs must be above the sensitivity limit of the observa- 
tions and must be clearly distinguishable from the Galac- 
tic background emission (Green 1991). Therefore, the data 
from the fainter parts of the loop (which have very low sur- 
face brightness) are not taken into account because it is very 
difficult to resolve them from the background. On the other 
hand, this would significantly reduce the brightnesses of en- 
tire loops and there is a general trend that fainter SNRs tend 
to be larger (Green 2005). Consequently, the S — Z) rela- 
tions that were applied could give overestimated values for 
the diameters. For the estimation of brightness temperatures 
over the spurs we had to take into account background ra- 
diation (Webster 1974). Borders enclosing the spurs are de- 
fined such as to clearly delineate a spur and its background. 
For the method of calculation see Borka (2007). 

The loops are believed to be SNRs whose material is 
expanding inside of bubbles of low density, created either 
by earlier SNR explosions or by strong stellar winds (Salter 
1983; McKee & Ostriker 1977 and references therein). 

The locations of the Galactic radio loops (and their spu- 
rs) in the sky, and the values of their average brightness tem- 
peratures and surface brightnesses, can be found in Borka 
(2007; Tables 1 and 2 therein). 

The temperatures at 1 GHz were calculated using the 
temperatures at 1420 MHz given by Borka (2007) and the 
spectral indices of Berkhuijsen (1973) for Loops I-IV, plus 
those calculated here for Loops V and VI. We calculated the 
surface brightness by applying the following relations: 



TiGHz/ri.42GHz - (L0/L42)-'3 

ElGHz = 2fcriGHz(lGHz/c)2 



(2) 
(3) 



where /? is the temperature spectral index. 

In order to determine distances for the Galactic radio 
loops I- VI, we applied the surface brightnesses from Borka 
(2007) to the recent E — _D relations of Urosevic (2002) and 
Arbutina et al. (2004). In addition, the values of the spec- 
tral indices of the main Galactic Loops I-IV were obtained 
directly from spectra of temperature versus frequency. To 
produce these spectra, with the 1420 MHz data from Borka 
(2007) (map of Reich & Reich 1986) we used data from 
Berkhuijsen (1973) at 820 MHz for Loops I-III, and at 404 
MHz for Loop IV. The survey at 820 MHz is by Berkhui- 
jsen (1972) and that at 404 MHz by Seeger et al. (1965). 
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Table 1 The diameters D (pc) and distances r (pc) to 
the main radio loops I-IV derived from the Yj — D relations 
given by Urosevic (2002) and Arbutina et al. (2004), respec- 
tively. 



Loop 


derived from (|4j 


derived from Q 


Z3(pc) 


r (pc) 


D(pc) 


r (pc) 


I 


169 ± 38 


99 ±25 


119 ± 17 


70 ± 12 


11 


187 ±42 


129 ± 33 


127 ± 18 


88 ± 16 


III 


167 ± 37 


150 ± 42 


118 ± 17 


106 ± 21 


IV 


310 ± 69 


462 ± 103 


172 ± 24 


257 ± 50 



The T — T method was used for the determination of spec- 
tral indices of Loops V and VI. 

Diameters were calculated using the E — D relations 
(Urosevic 2002; Arbutina et al. 2004): 

SiGHz = 8.43 • 10-^^ D-^-^^ (4) 
EiGHz = 2.5 • 10-15^^-3.4 

(EiGHz in Wm-^ Hz^ Si^^ and D in pc) and from these 
the distances to the centers of the shells were calculated us- 
ing the relation: 

r = D/{2sm0) (6) 

with angular radii (9) taken from Milogradov-Turin & Uro- 
sevic (1997). In Urosevic (2002) the Galactic loops are ad- 
ded to the set of 157 calibrators, and relation (|4|i was de- 
rived. In Arbutina et al. (2004) a good S - L» relation (01 
was derived for the M82 starburst galaxy and they suggest 
that a similar relation might exist for Galactic SNRs associ- 
ated with giant molecular clouds. 
In addition, we have the relation: 

Tb = Ki^-I^ (7) 

where (3 is the spectral index and K isa constant. Knowing 
the values of brightness temperature at two frequencies, we 
were able to calculate spectral indices. 

For Loops V and VI, spectral indices were calculated 
using T — T plots. These were made between temperatures 
at 1420 MHz (Reich & Reich 1986) and 408 MHz (Haslam 
et al. 1982). 

By the relation: 

drb(i^i)/dTb(i^2) - {1^1/1^2)-^ (8) 

we can determine the temperature spectral index, (3, from 
the slope of straight-line fitted to a T — T plot. Using the 
above relations, diameters and distances were calculated for 
Loops V and VI. 

4 Results 

The results for the four main loops are presented in Tables 
[T]and|2] The diameters and distances to the radio loops de- 
rived from the Y, — D relations given by Urosevic (2002) 



Table 2 The spectral indices of the main radio loops. 



Loop 


/3 


I 


2.69 ± 0.05 


II 


2.91 ± 0.05 


III 


2.59 ± 0.04 


IV 


2.65 ± 0.09 



Table 3 The diameters D (pc) and distances r (pc) to the 
radio loops V and VI derived from the Y — D relations given 
by Urosevic (2002) and Arbutina et al. (2004), respectively. 



Loop 


derived from (O 


derived from ([5j 


D(pc) 


r (pc) 


D(pc) 


r (pc) 


V 


235 ± 36 


124 ± 22 


146 ± 22 


77 ± 14 


VI 


255 ± 39 


134 ± 23 


153 ± 23 


80 ± 14 



Table 4 The spectral indices of radio loops V and VI, 
between 1420 and 408 MHz. 



Loop 


Spur 






V 


Taurus 


2.89 ±0.10 


2.90 ± 0.11 


Pisces 


2.63 ±0.12 


Pegasus 


3.12 ±0.09 


VI 


Leo 


2.91 ±0.13 


2.77 ±0.13 


Cancer 


2.71 ±0.12 



and Arbutina et al. (2004) are shown in Table [T] The spec- 
tral indices for the four main radio loops, calculated as de- 
scribed above, are shown in Table |2l Spectra were gener- 
ated using mean temperatures at two different frequencies 
as needed for making the simplest linear fits. The two fre- 
quencies are 820 and 1420 MHz for Loops I-III, and 404 
and 1420 MHz for Loop IV. All earlier determinations of 
the radio loop spectral indices, except for this paper and 
Borka (2007), were based on T — T methods. Besides, the 
values obtained (all greater than 2.2) confirm a non-thermal 
origin for the emission from the radio loops. The values for 
the brightness temperature spectral indices of the loops are 
rather steep (about 2.5 < /3 < 3.0). This is at the high end of 
the spectral index distribution for SNRs (Clark & Caswell 
1976). 

The diameters of, and distances to. Loops V and VI de- 
rived using the Y — D relations given by Urosevic (2002) 
and Arbutina et al. (2004) are presented in Table [3] Repre- 
sentative T — T plots between 1420 and 408 MHz are given 
in Figures |2] and [3] with one plot for each loop. Using simi- 
lar plots, we determined spectral indices for all spurs which 
create Loops V and VI. The spectral indices of Loops V 
and VI are shown in Table |4] The values of spectral index 
derived from T — T plots for these two loops are obtained 
independently of the values for the indices derived from the 
mean brightness temperatures in Borka (2007). 

We note that the values of the spectral indices in Tables|2] 
and|4]are slightly different to those derived by Borka (2007). 
However, the results agree within the absolute error limits. 
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T(K) 1420 MHz 

Fig. 2 The r - r plot for the Taurus spur (Loop V) at 
b = -13°. T (K) at 408 MHz is plotted versus T (K) at 
1420 MHz, for the longitude range from 178° to 189°. 
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T(K) 1420 MHz 

Fig. 3 The T - T plot for the Cancer spur (Loop VI) at 
6 = 9°. T (K) at 408 MHz is plotted versus T (K) at 1420 
MHz, for the longitude range from 196° to 207°. 

5 Discussion 

We will assume below that all the radio loops discussed here 
are SNRs (Berkhuijsen, Haslam & Salter 1970; Berkhuijsen 
1971; Shklovskii & Sheffer 1971; Salter 1983). The geome- 
tries of Loops V and VI are presented in Figure [1] It is seen 
that Loop V is projected entirely within Loop VI indicating 
a possible close connection between these loops. Here, we 
will consider three possibilities for such a connection: 

(i) Milogradov-Turin & Urosevic (1997) considered the po- 
ssibility that all the spurs of Loops V and VI define a unique 
small circle corresponding to a single feature - Loop Vh-VI. 
However, their geometric fit did not provide reliable evi- 
dence for this hypothesis. 

(ii) a second possibility is that Loop V represents a "wrin- 
kle" on the shell of Loop VI, tracing part of a small circle 
on that shell. If so, this small circle is oriented such that at 



one point (i.e. I = 150°, b = -44°) it meets the (differ- 
ent) small circle representing the locus of tangential points 
from the Sun to the Loop VI shell. Borka (2007) performed 
the first spectral investigation of Loops V and VI using the 
mean brightness temperatures at three frequencies (1420, 
820 and 408 MHz) to form the linear spectra for these loops. 
As shown in Table 1 of Borka (2007), spurs of these loops 
belong to different and not overlapping temperature inter- 
vals at all three frequencies. This suggests that these spurs 
form two independent features, namely Loops V and VI. 
(iii) a similar situation to that detailed in (i) and (ii) but con- 
cerning Loops I and IV, was considered earlier by a number 
of authors (Large et al. 1962; Davies 1964; Guidice 1971). 
However, Borken and Iwan (1977) and Iwan (1980) pro- 
posed that the Loop I SNR is very old, well into its radia- 
tive phase, and is being reheated by the intersecting shock 
of Loop IV. They concluded that Loop IV has a low den- 
sity shell, is not yet in the radiative cooling phase, and has 
not lost much energy before reheating the Loop I cavity. By 
analogy, and given the similar distances and diameters de- 
rived for these two loops (TableO, we assume that the Loop 
V SNR exploded within the previous (low density) cavity 
left by the Loop VI SNR. In view of the similarity of the 
paths of the Loops V and VI, one possible scenario is sug- 
gested by the recent paper by Wolleben (2007). This author 
proposes that the emission of two independent synchrotron- 
emitting shells make up the feature known as Loop I (the 
North Polar Spur). While the presence of more than one 
supernova remnant has been invoked previously to explain 
some of the features of Loop I, (and especially its strong 
X-ray emission), WoUeben's model invokes the presence of 
a second, very old remnant to explain certain features he 
finds in the large-scale distribution of polarized continuum 
intensity. 

We have used recent Y, — D relations to determine the 
distances of the various loops. Our results demonstrate that 
the distances from Arbutina et al. (2004) do not agree with 
those derived independently by Berkhuijsen (1986). This is 
probably because relation (|5]l is not applicable to the Loops, 
as these are thought to be situated in a low density environ- 
ment. An analysis as to how the Y, — D relation depends 
on the density of the interstellar medium is given in Ar- 
butina & Urosevic (2005). Therefore, we conclude that the 
radio loops are evolving in a low density environment. How- 
ever, we emphasize once again our assumption of the SNR 
explanation for the radio loops. The distances and diame- 
ters obtained using the Y, ~ D relation of Urosevic (2002) 
show good agreement when compared with those in Borka 
(2007), which were calculated from the Y — D relation de- 
rived by Case & Bhattacharya (1998). However, the relation 
of Urosevic is preferred to that of Case & Bhattacharya be- 
cause it adds the Galactic loops to the set of 157 calibra- 
tors. Also, the relation of Case & Bhattacharya shows better 
agreement with the diameters of SNRs at known distances 
(Berkhuijsen 1986) than does that of Arbutina et al. (2004). 
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Fig. 4 The surface brightness to diameter diagram from 
Berezhko & Volk (2004), with values for Loops I- VI added. 
Three different densities for the ISM (Nh = 3, 0.3 and 0.003 
cm^'^) are presented, plus three values for the SN explosion 
energy (£:sn = 0.25, 1 and 3 x lO^i erg). 

Adopting diameters for the main loops from Berkhui- 
jsen (1986) and the diameters for Loops V and VI as es- 
timated in this paper from relation (|4|i, we can use the di- 
agram of Berezhko & Volk (2004) to determine the ini- 
tial explosion energy that would correspond to these loops, 
plus the density of the ambient interstellar medium in their 
neighborhood. The — D diagram at 1 GHz taken from 
Berezhko & Volk (2004) with the derived values for the six 
loops superposed, is shown in Figure 2] We identify Loops 
I-IV by circles and Loops V and VI by stars. For Loops V 
and VI, this is a preliminary illustration. For additional evi- 
dence of their SNR origin it would be necessary to calculate 
distances (diameters) independently of the E — D relation. 
Then, from this E — diagram, it would be possible to in- 
vestigate the compatibility of the properties of Loops V and 
VI with their being old SNRs. From Figure|4] it can be con- 
cluded that the surrounding density of the radio loops is low 
and that the explosion energy is high. Further, the locations 
of the Loops on this Y, — D diagram agree with the general 
distribution and this strengthens the idea that they are old 
SNRs. 

The main unsolved problems for the SNR hypothesis of 
the radio loops are: (i) the high X-ray flux from the interior 
of Loop I, and (ii) the high degree of spherical symmetry 



of all the loops. Both issues are related to the evolutionary 
status of radio loops. If the loops are SNRs, they are old 
and dimensionally immense (see Tables [T] and O. High X- 
ray emission and spherical symmetry are not expected from 
such SNRs. 

The X-ray contours follow the North Polar Spur (NPS), 
but they are located in the interior region of Loop I (Sofue 
2000). We believe that Loop I was in a late adiabatic phase, 
not so long ago, and in that phase X-rays could be gener- 
ated. At present, the NPS is probably in an early radiative 
phase of expansion in which a thick low temperature shell 
is generated (Chevalier 1974). The X-ray gas is thermody- 
namically stable (the hot phase of the ISM) and can exist for 
a million years, implying that X-rays can be observed from 
the interior of Loop I. In a thicker uniform medium, a slow 
shock can conserve its own approximately spherical form, 
but only for a very short time interval. Due to the formation 
of a thick shell during expansion in the radiative phase (Cox 
& Reynolds 1987), the velocity of the shock wave rapidly 
declines and the shell loses its own identity. Generally, the 
radiative phase of evolution is short in comparison to the 
adiabatic phase (e.g. McKee & Ostriker 1977). 

If a SNR is expanding in a hot bubble made by strong 
winds or earlier supernova explosions, it has a chance to 
conserve an approximately spherical form because of the 
uniformly distributed matter formed by the evaporation of 
denser clouds (McKee & Ostriker 1977). Hanayama & To- 
misaka (2006) concluded that a spherically symmetric sha- 
pe can be conserved if an old SNR evolves in an environ- 
ment with a weak magnetic field. In this case, the evolution 
is similar to that of a non-magnetized SNR. In the pressure- 
driven expansion stage, a SNR evolving in a medium with 
a magnetic field of 1 /iG keeps its spherically symmetric 
shape (Hanayama & Tomisaka 2006). 

The above discussion suggests that the radio loops may 
have a SNR origin, but a definite answer needs much addi- 
tional work. 

6 Summary 

The main results of this paper may be summarized as fol- 
lows: 

(i) We have estimated the distances of the main Galactic ra- 
dio loops and Loops V and VI from the northern-sky radio 
continuum survey at 1420 MHz (Reich & Reich 1986) us- 
ing two different E — _D relations. 

(ii) The estimated distances of the main radio loops derived 
using the 1420 MHz data are in good agreement with the 
earlier results (e.g. Berkhuijsen 1986). 

(iii) From the spectral index analysis we can confirm that 
the emission from the radio loops is non-thermal in origin. 

(iv) We have demonstrated the probable existence of Loops 
V and VI. The distances and spectral indices of Loop V and 
Loop VI are estimated. In general, they display shell struc- 
ture and have non-thermal spectra, similar to the main loops 
I-IV which are usually assumed to be SNRs. 
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